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Abstract: In the framework of three-point QCD sum rules, the form factors for the semileptonic 
decays of — > B s (B*)l + Vi are calculated with account for the Coulomb-like a s /i;-corrections in the 
heavy quarkonium. The generalized relations due to the spin symmetry of HQET/NRQCD for the 
form factors are derived at the recoil momentum close to zero. The nonleptonic decays are studied 
using the assumption on the factorization. The B c meson lifetime is estimated by summing up the 
dominating exclusive modes in the c — > s transition combining the current calculations with the 
previous analysis of b — > c decays in the sum rules of QCD and NRQCD. 
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1. Introduction 

For better understanding and precise measuring 
the weak-action properties of heavy quarks, gov- 
erned by the QCD forces, we need as wide as pos- 
sible collection of snapshots with hadrons, con- 
taining the heavy quarks. Then we can provide 
the study of heavy quarks dynamics by testing 
the various conditions, determining the forming 
of bound states as well as the entering of strong 
interactions into the weak processes. So, a new 
lab for such investigations is a doubly heavy long- 
lived quarkonium B c recently observed by the 
CDF Collaboration [Q for the first time. 

This meson is similar to the charmonium and 
bottomonium in the spectroscopy, since it is com- 
posed by two nonrclativistic heavy quarks, so 
that the NRQCD approach [H is well justified 
to the system. The modern predictions for the 
mass spectra of be levels were obtained in refs. 
H in the framework of potential models and lat- 
tice simulations. The measured value of B c mass 
yet has a large uncertainty Mb c = 6.40 ± 0.39 ± 
0.13 GeV, in agreement with the theoretical ex- 
pectations. 

The measured B c lifetime 

t[B c ] = 0.46±°i| ± 0.03 ps, 
agrees with the estimates obtained in the frame- 



work of both the OPE combined with the eval- 
uation of hadronic matrix elements in NRQCD 
|4|, ||, [| and potential quark models, where one 
has to sum up the dominating exclusive modes to 
calculate the total B c width 0, ||, iope,pm[-8 c ] = 
0.55 ± 0.15 ps. The accurate measurement of B c 
lifetime could allow one to distinguish various pa- 
rameter dependencies such as the optimal heavy 
quark masses, which basically determine the the- 
oretical uncertainties in OPE. 

At present, the calculations of B c decays in 
the framework of QCD sum rules were performed 
in §, [uj 0, ||. The authors of |, [RJ got 
the results, where the form factors are about 3 
times less than the values expected in the po- 
tential quark models, and the semileptonic and 
hadronic widths of B c are one order of magnitude 
less than those in OPE. The reason for such the 
disagreement was pointed out in JO] and stud- 
ied in [jL2[ : in the QCD sum rules for the heavy 
quarkonia the Coulomb-like corrections are sig- 
nificant, since they correspond to summing up 
the ladder diagrams, where a s /v is not a small 
parameter, as the heavy quarks move nonrela- 
tivistically, v <C 1. The Coulomb rescaling of 
quark-quarkonium vertex enhances the estimates 
of form factors in the QCD sum rules for the 
£?+ — > il>(rj c )l + v decays. In the soft limit 
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v\ ■ v 2 — > 1 at vi V2, where Ux,2 denote the 
four-velocities of initial and recoil mesons, was 
considered, and the generalized spin symmetry 
relations were obtained for the B c — > ip(rj c ) tran- 
sitions: four equations, including that of 
Moreover, the gluon condensate term was cal- 
culated in both QCD and NRQCD, so that it 
enforced a convergency of the method. 

In the present paper we calculate the B c de- 
cays due to the c — ► s weak transition in the 
framework of QCD sum rules, taking into ac- 
count the Coulomb- like a s /u-corrections for the 
heavy quarkonium in the initial state. In the 
semileptonic decays the hadronic final state is 
saturated by the pseudoscalar B s and vector B* 
mesons, so that we need the values of their lep- 
tonic constants entering the sum rules and de- 
termining the normalization of form factors. For 
this purpose, we reanalyze the two-point sum 
rules for the B mesons to take into account the 
product of quark and gluon condensates in addi- 
tion to the previous consideration of terms with 
the quark and mixed condensates. We demon- 
strate the significant role of the product term for 
the convergency of method and reevaluate the 
constants fs as well as fs, ■ Taking into account 
the dependence on the threshold energy E c of 
hadronic continuum in the bs system in both the 
value of /b s extracted from the two-point sum 
rules and the form factors in the three-point sum 
rules, we observe the stability of form factors ver- 
sus E c , which indicates the convergency of sum 
rules. 

The spin symmetries of leading terms in the 



lagrangians of HQET 14 for the singly heavy 
hadrons (here B { s * ] ) and NRQCD § for the dou- 
bly heavy mesons (here B c ) result in the rela- 
tions between the form factors of semileptonic 
B c — > Bg decays. We derive two generalized 
relations in the soft limit Vi • v% 1: one equa- 
tion in addition to what was found previously in 
rcf.|jl|. The relations are in a good agreement 
with the sum rules calculations up to the accu- 
racy better than 10%, that shows a low contri- 
bution of next-to- leading l/mQ-tcrms. 

We perform the numerical estimates of semi- 
leptonic B c widths and use the factorization ap- 
proach jlB) to evaluate the nonleptonic modes. 
Summing up the dominating exclusive modes, we 



calculate the lifetime of B c , which agree with the 
experimental data and the predictions of OPE 
and quark models. We discuss the preferable pre- 
scription for the normalization point of nonlep- 
tonic weak lagrangian for the charmed quark and 
present our optimal estimate of total B c width. 
We stress that in the QCD sum rules to the given 
order in a s , the uncertainty in the values of heavy 
quark masses is much less than in OPE. This fact 
leads to a more definite prediction on the B c life- 
time. 

2. Three-point sum rules 

The hadronic matrix elements for the semilep- 
tonic B c (pi) — > B s {p2) decays can be written 
down as follows: 

(B S \V^\B C ) = f+( Pl + p 2 )^ + /_<,„, (2.1) 

±{B*MB C ) = iF v e^ aP e*»( Pl + p») a q fi , 

\{B:\A,\B C ) = F A e; + F£(e*-p 1 )(p 1 + p 2 ), 

where = {pi -p%)n and e M = e^(p 2 ) is the 
polarization vector of B* meson. and are 
the flavour changing vector and axial electroweak 
currents. Following the standard procedure for 
the evaluation of form factors in the framework 
of QCD sum rules jl6), we consider the three- 
point functions, say, 

U^ Pl ,p 2 ,q 2 ) = l 2 J dxdye^-*-^ ■ 
(0\T{q 1 (x) j5 q 2 (x),V^Q),b(yh5c(y)}\0), 

where qi(x)^q 2 {x) and q~i(x)j l/ q 2 (x) denote in- 
terpolating currents for B s and B*, correspond- 
ingly. 

The Lorentz structures in the correlators can 
be written down as 11^ = 11+ (pi + p 2 )^ + Il-g^. 
The form factors f± are determined from the am- 
plitudes II±, respectively. 

The leading QCD term is a triangle quark- 
loop diagram, for which we can write down the 
double dispersion representation at q 2 < 

U pert (n 2 r? n 2 ) - 1 . 

dsids 2 + subtractions, 



P 



pert 



si,s 2 ,q 2 ) 



{si - p\){s 2 - p 2 ) 
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where the limits of integration region and the 
spectral densities are given in Q . 

The physical spectral functions are generally 
saturated by the ground hadronic states and a 
continuum starting at some effective thresholds. 

3 




Figure 1: The ladder diagram of the Coulomb- like 
interaction. 

For the heavy quarkonium be, where the rel- 
ative velocity of quark movement is small, an es- 
sential role is taken by the Coulomb- like a s /v- 
corrections. They are caused by the ladder dia- 
gram, shown in Fig. [l]. This leads to the finite 
renormalization for pi p2[, so that p\ — Cp%, 



C = 



*£(o)| 



I* 



free 
be 



(0)1 



47ra s 47ra s 
(l-exp{ — I)- 1 , 



3r 



3r 



where v is the relative velocity of quarks in the 
6c-system, v = Jl- , 4 , m "" ^^. 

3. Numerical estimates 

We evaluate the form factors in the scheme of 
spectral density moments. This scheme is not 
strongly sensitive to the value of the 6c-system 
threshold energy, and we put E b c c = 1.2 GeV. 
The two-point sum rules for the B c meson with 
account for the Coulomb-like corrections give 

a c s (bc) = 0.45, 

which corresponds to / Bc =400 MeV 0. The 
quark masses are fixed by the calculations of lep- 
tonic constants and /t in the same order 
over a s . The requirement of stability in the sum 
rules including the contributions of higher ex- 
citations, results in quite an accurate determi- 
nation of masses m c = 1.40 ± 0.03 GeV and 
nib = 4.60 ±0.02 GeV, which are in a good agree- 
ment with the recent estimates in 118L where 



the quark masses free off a renormalon ambiguity 
were introduced. 

The leptonic constant for the B s meson is ex- 
tracted from the two-point sum rules. The Borel 
improved sum rules for the B meson leptonic con- 
stant p9| have the following form: 



u>o(f0 



f B M B e- K ^ T =K 2 ^C{p) / duju?e 



,2 J2 



16 



2 A. 



S G 2 », 



288 x 7T 



where the K-factor is due to a s -corrections |19| . 
We find that NLO corrections to the leptonic 
constant are about 40%. Using the Pade approx- 
imation, we find that higher orders corrections 
can be about 30%. So, we hold the K factor in 
conservative limits 1.4-5-1.7. ft is quite reason- 
able to suppose its cancellation in evaluating the 
semileptonic form factors due to the renormaliza- 
tion of heavy- light vertex in the triangle diagram. 
In the limit of semi-local duality [^0[ |2l| t — > 
we get the relation: A(/i) = | w (/i). We in- 
troduce the renormalization invariant quantities 

For 



<7ual = C-^(p) WoO*), A™ = | W 



Q.dual 



A r / n , we have A r f n , 

dual dual 



= M B -m b = 0.63 GeV, 
and we obtain that in the semi-local duality the 



threshold energy 



O.dual 



= 0.84 GeV. Neglect- 



ing the quark condensate term in the leptonic 
constant we have f%M B = K 2 JrK ; X oi ) 3 . In 
the general Borel scheme for fg we have to con- 
sider the stability at r ^ with the extended re- 
gion of resonance contribution. We expect, that 
the sum rules with the redefined ui ren and A rera 
have a stability point at r ~ -r. The results are 
in a good agreement with the semi-local dual- 
ity if the threshold energy of continuum equals 
E c = 1.1 -T- 1.3 GeV (see Fig. |2|, where the over- 
all K-factor was ignored). So, we find the E^J 2 - 
dependence of fg V M B , whereas the contribu- 
tion of condensate is numerically suppressed, as 
expected from the semi-local duality. Multiply- 
ing the result taken from Fig. ||, by the K-factor 
we find the value fg — 140 -j- 170 MeV, which is 
in a good agreement with the recent lattice re- 
sults p2) and the estimates in the QCD SR by 



other authors |23 



put 



For the vector B* meson constant fg* we 
1.11 (see @, H). For the leptonic 



Ib 
Jb 
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f B VM, GeV 3/2 A, GeV 
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t, GeV' 1 t, GeV" 1 

Figure 2: The leptonic constant of B meson and the b-quark binding energy in the semi- local duality sum 
rules (dashed curve) and in the general Borel scheme (solid line) with the corrected value of A, which improves 
the stability of result obtained in the semi-local duality. 



constant of B s meson we explore the following 
relation = 1.16, which expresses the flavor 
SU(3)-symmetry violation for B mesons pjj . 

We have investigated the dependence of form 
factors on the bs threshold energy of continuum 
in the range E c = 1.1 -j- 1.3 GeV, so that the op- 
timal choice for the bs system threshold energy is 
1.2 GeV. In Table [l] we present the results of sum 
rules for the form factors. Comparing with the 
estimates in the framework of potential models 
||, |2f| , we find a good agreement of estimates in 
the QCD sum rules with the values in the quark 
model. 



/+ 


/- 


iV,GeV _1 


F A ,GeV 


1.3 


-5.8 


1.1 


8.1 



Table 1: The form factors of B c decay modes into 
the B B and mesons at q 2 — 0. 

The accuracy of sum rules under considera- 
tion is basically determined by the variation of 
heavy quark masses. Indeed, the significant a s 
correction to the leptonic constant of B s meson 
should cancel the same factor for the renormal- 
ization of quark- meson vertex in the triangle dia- 
gram. The dependence on the choice of threshold 
energy in the 6s-channel can be optimized and, 
hence, minimized. The variation of threshold en- 
ergy in the fee-channel give the error less than 1%. 
The effective coulomb constant is fixed from the 
two-point sum rules for the heavy quarkonium, 
and its variation is less than 2%, which gives the 
same uncertainty for the form factors. The heavy 



quark masses are determined by the two-point 
sum rules for the heavy quarkonia, too. How- 
ever, their variations result in the most essential 
uncertainty. Summing up all of mentioned vari- 
ations we estimate 6f/f~ 5%. 

The semileptonic widths calculated in the 
QCD sum rules are presented in Table 0. 



mode 


r, 1CT 14 GeV 


BR, % 


B s e + v e 


5.8 


4.0 


B* s e+y e 


7.2 


5.0 



Table 2: The widths of semileptonic B c decay modes 
and the branching fractions calculated at tb c = 0.46 
ps. 



4. The symmetry relations 

At the recoil momentum close to zero, the heavy 
quarks in both the initial and final states have 
small relative velocities inside the hadrons, so 
that the dynamics of heavy quarks is essentially 
nonrelativistic. This allows us to use the com- 
bined NRQCD /HQET approximation in the stu- 
dy of mesonic form factors. The expansion in the 
small relative velocities leads to various relations 
between the form factors due to the spin symme- 
try of effective lagragians to the leading order. 
Solving these relations results in the introduc- 
tion of an universal form factor (an analogue of 
the Isgur-Wise function) at q 2 — > q^ax- 
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We have derived the symmetry relations for 
the following form factors: 

/+(cf .M 2 -4Mx)- 

/_(cf -M 2 + c% -Mi) = 0, 

■ c v - 2c, ■ F v MxM 2 = 0, (4.1) 

F A c^ + c e -Mi (/+ + /_) = 0, 

where M\ = m c + nib, M-i = ~m s + mi, and 



c e = -2, 
cy = -1-5- 



l 



= 1 - B 



2m c ' 



4 



l+B 



2m c 
m b 
2m c 



(4.2) 



Equating the second relation in (4.1), for exam- 
ple, we obtain 



B = — 



2m c + mb 4mb(m c + rrib)Fv 



2m c 



F A 



10.0, 



where all form factors are taken at q^ax- Sub- 
stituting B in first and third relations, we get 
/+ ~ 2.0 and /_ ~ —8.3. These values have to 
be compared with the corresponding form factors 
obtained in the QCD sum rules: f + (q 2 nax ) = 1.8 
and f-(qmax) = —8-1, where we suppose the pole 
like behaviour of form factors. Thus, we find that 



in the QCD sum rules, relations (4.1) are valid 
with the accuracy better than 10% at q 2 = q^ax- 
The deviation could increase at q 2 < q^ax be- 
cause of variations in the pole masses governing 
the evolution of form factors. However, in B^ — ► 
Bg 'l + v decays the phase space is restricted, so 
that the changes of form factors are about 50%, 
while their ratios develop more slowly. 



5. Nonleptonic decays and the life- 
time 

The hadronic decay widths can be obtained on 
the basis of assumption on the factorization for 
the weak transition between the quarkonia and 
the final two-body hadronic states. For the non- 
leptonic decay modes the effective Hamiltonian 
can be written down as 



H, 



eff 



G F 
2^/2 



v cs v: d {c+(ix)o + + c^)0-}, 



where 0± = {uij„(l - ^)di)(sjY{^ ~ 7s)cj) ± 
(iiij u (l - j 5 )dj)(s~iY(l - lh)Cj), and the factors 
C±(a*) account for the strong corrections to the 
corresponding four-fermion operators caused by 
hard gluons. The review on the evaluation of 
C±(aO can be found in pf|. The results are col- 
lected in Table [|. 



mode 


T, 10- 14 GeV 


BR, % 


B s n+ 


15.8 a\ 


17.5 


B sP + 


6.7 a\ 


7.4 


B*TT+ 


6.2 a\ 


6.9 


B* sP + 


20.0 a\ 


22.2 



Table 3: The widths of dominant nonleptonic B c de- 
cay modes due to c — » s transition and the branch- 
ing fractions calculated at tb c = 0.46 ps. We put 
ai=1.26. 

In the parton approximation we could ex- 



pect F[B+ -> B. 



(*) 



-light hadrons] = (2C? (fj,) 



B^e+v P 



which results in the 



estimate very close to the value obtained as the 
sum of exclusive modes at fi > 0.9 GeV. The de- 
viation between these two estimates slightly in- 
crease at ^ < n < 0.9 GeV. Concerning the 
comparison of hadronic width summing up the 
exclusive decay modes with the estimate based 
on the quark-hadron duality, we insist that the 
deviation between these two estimates is unessen- 
tial since it is less that 10%. 

We estimate the lifetime using the fact that 
the dominant modes of the B c meson decays are 
the c — > s, b — > c transitions with the Bs and 
J/t/j, r)c final states respectively, and the elec- 
troweak annihilation . 

The method for the calculation of multi-par- 
ticle branching fractions was offered by Bjorken 
in his pioneering paper devoted to the decays 
of hadrons containing heavy quarks |27j]. In or- 
der to estimate the contribution of non-resonant 
3ir modes of B c decays into B^ we use this 
technique, i.e. the Poisson distribution with the 
average value corrected to agree with the non- 
resonant 37r-modes in the decays of D mesons. 
We have found BR(B+ -> si* ) (37r)+) « 0.2%, 



The b — > ccs transition is negligibly small in the B c 
decays because of destructive Pauli interference for the 
charmed quark in the initial state and the product of de- 
cay §j. 
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while BR(B+ Sr ) (2^)+| non _ rC sona„t) « 3%. 
We see that the neglected modes contribute to 
the total width of B c as a small fraction in the 
limits of uncertainty envolved. 



t, ps 




0.8 1 1.2 1.4 1.6 

fi, GeV 



Figure 3: The dependence of B c meson lifetime on 
the scale /i in the effective Hamiltonian (^). The 
shaded region shows the uncertainty of estimates, the 
dark shaded region is the preferable choice as given 
by the lifetimes of charmed mesons. The dots repre- 
sent the values in the OPE approach. 

The width of beauty decay in the sum rules 
was derived using the similar methods in fil| : 
T{B+ cc + X) = (28 ± 5) ■ 1(T 14 GeV. The 
width of the electroweak annihilation is taken 
from || as 12 • 1CT 14 GeV. 

In Fig. [| we present the B c meson lifetime 
calculated in the QCD SR under consideration. 
We also show the results of the lifetime evalua- 
tion in the framework of Operator Product Ex- 
pansion in NRQCD §, §. 

In contrast to OPE, where the basic uncer- 
tainty is given by the variation of heavy quark 
masses, these parameters are fixed by the two- 
point sum rules for bottomonia and charmonia, 
so that the accuracy of SR calculations for the 
total width of B c is determined by the choice of 
scale [i for the hadronic weak lagrangian in de- 
cays of charmed quark. We show this dependence 
in Fig. ||, where ^ < /i < m c and the dark 
shaded region corresponds to the scales preferred 
by data on the charmed meson lifetimes. The dis- 
cussion on the optimal choice of scale in hadronic 
decays is addressed in Q . We suppose that the 
preferable choice of scale in the c — > s decays of 
B c is equal to fi B — fi c i ■ fi cS w (0.85 GeV) 2 , and 



at ai(pB c ) = 1-20 in the charmed quark decays 
we predict t[B c ] = 0.48 ± 0.05 ps. 

6. Conclusion 

We have investigated the semileptonic decays of 
B c meson due to the weak decays of charmed 
quark in the framework of three-point sum rules 
in QCD. We have pointed out the important role 
played by the Coulomb- like a s /w-corrections. As 
in the case of two-point sum rules, the form fac- 
tors are about three times enhanced due to the 
Coulomb renormalization of quark-meson vertex 
for the heavy quarkonium B c . We have studied 
the dependence of form factors on the thresh- 
old energy, which determines the continuum re- 
gion of bs system. The obtained dependence has 
the stability region, serving as the test of con- 
vergency for the sum rule method. The HQET 
two-point sum rules for the leptonic constant fg s 
and Jb- have been reanalyzed to introduce the 
term caused by the product of quark and gluon 
condensates. This contribution essentially im- 
proves the stability of SR results for the leptonic 
constants of B mesons, yielding: fg = 140 -f- 170 
MeV. 

We have studied the soft limit for the form 
factors in combined HQET/NRQCD technique 
at the recoil momentum close to zero, which al- 
lows us to derive the generalized relations due to 
the spin symmetry of effective lagrangian. The 
relations are in a good agreement with the full 
QCD results, which means that the corrections 
to the form factors in both relative velocity of 
heavy quarks inside the bc quarkonium and the 
inverse heavy quark masses are small within the 
accuracy of the method. 

Next, we have studied the nonleptonic de- 
cays, using the assumption on the factorization 
of the weak transition. The results on the widths 
and branching fractions for various decay modes 
of B c are collected in Tables. 

Finally, we have estimated the B c meson life- 
time, and showed the dependence on the scale 
for the hadronic weak lagrangian in decays of 
charmed quark t[B c ] = 0.48 ± 0.05 ps. Our es- 
timates are in a good agreement with the the- 
oretical predictions for the lifetime in both the 
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potential models and OPE as well as with the 
experimental data. 
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of RFBR 99-02-16558 and 00-15-96645. 
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